Sporadic randomness, Maxwell's Demon and the Poincare recurrence times 
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In the case of fully chaotic systems the distribution of the Poincare recurrence times is an exponential 
whose decay rate is the Kolmogorov-Sinai (KS) entropy. We address the discussion of the same 
problem, the connection between dynamics and thermodynamics, in the case of sporadic randomness, 
using the Manneville map as a prototype of this class of processes. We explore the possibility of 
relating the distribution of Poincare recurrence times to "thermodynamics" , in the sense of the KS 
entropy, also in the case of an inverse power law. This is the dynamic property that Zaslavsly ( 
G.M. Zaslavsky, Physics Today 52(8), 39 (1999)) finds to be responsible for a striking deviation 
from ordinary statistical mechanics under the form of Maxwell's Demon effect. 

We show that this way of establishing a connection between thermodynamics and dynamics is valid 
only in the case of strong chaos, where both the sensitivity to initial conditions and the distribution 
of the Poincare recurrence times are exponential. In the case of sporadic randomness, resulting at 
long times in the Levy diffusion processes, the sensitivity to initial conditions is initially a power 
law, but it becomes exponential again in the long-time scale, whereas the distribution of Poincare 
recurrence times keeps, or gets, its inverse power law nature forever, including the long-time scale 
where the sensitivity to initial condition becomes exponential. 

We show that a non-extensive version of thermodynamics would imply the Maxwell's Demon effect 
to be determined by memory, and thus to be temporary, in conflict with the dynamic approach to 
Levy statistics. The adoption of heuristic arguments indicate that this effect is possible, as a form 
of genuine equilibrium, after completion of the process of memory erasure. 
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I. INTRODUCTION 

In a lucid discussion Lebowitz |lj has recently restated 
the point of view of Boltzmann to establish the micro- 
scopic origin of irreversible macroscopic behavior. In his 
view the adoption of the laws of big numbers is essen- 
tial, and the role of deterministic chaos becomes impor- 
tant only if it applies to a macroscopic number of non- 
interacting particles. According to Lebowitz mixing and 
crgodicity are notions "unnecessary, misguided and mis- 
leading" . In other words, this opinion reflects the convic- 
tion, mirrored by the handbooks of statistical mechanics, 
that the unification of mechanics and thermodynamics 
rests on iV — > oo, where N denotes the number of de- 
grees of freedom of the system under study. 

This point of view has to be contrasted with a new 
view about the origin of thermodynamics [g| produced 
by the increasing interest in nonlinear dynamics. Within 
this new perspective the sensitivity to the initial condi- 
tions, rather than the large number of degrees of freedom 
involved, is playing the crucial role of generator of ther- 
modynamics. Zaslavsky [|| has recently discussed the 
phenomenon of the Poincare recurrences without using 
the traditional argument that these recurrences occur in 
a so extended time scale as to lie beyond the observa- 
tion range. Using the Bernouilli shift map, as a dynamic 



prototype of strong chaos, Zaslawsky proves [[| that the 
distribution of Poincare recurrence times, Pn(t), obeys 
the following prescription: 



P R (t) oc exp(-h KS t), 



(1) 



where Hks denotes the Kolmogorov-Sinai (KS) entropy 
[p| . Zaslavsky makes the interesting remark that Eq. ( l]) 
establishes a nice connection between mechanics, the left 
hand side of this equation, and thermodynamics, the 
right hand side of the same equation. This is so be- 
cause has has an entropic significance Q . It is also well 
known that the KS entropy is related to the Lyapunov 
coefficient and thus to the sensitivity to initial conditions, 
as a consequence of the Pesin theorem [Q. To make the 
connection between Pr(£) and the Lyapunov coefficient 
perspective more transparent, it is convenient to intro- 
duce also the property: 



£(t) = lim 



Ax{t) 



Axg^O Ax(0) 



(2) 



where Ax(t) denotes the difference between the space 
coordinates of two trajectories, moving from two distinct 
initial conditions at a distance |Aa;(0)| the one from the 
other. It is evident that after identifying has with the 
Lyapunov coefficient A, we can rewrite Eq.(|l|) in the form 



1 



where 



P R (t) OC £(-t) = — y 



= exp(Xt). 



(3) 



(4) 



At a first sight, this point of view seems to conflict 
with that of Lebowitz. In this paper we shall argue, with 
a heuristic approach, that it is not so, and that a sort of 
compromise between the two views can be found. Here 
we want to limit ourseleves to noticing that the observa- 
tion frequently made that the occurrence of strong chaos 
is exceptional, does not necessarily mean that thermody- 
namics can only rest on the (N — ► oo)-perspective. It 
might only imply that the process of transition to equi- 
librium and the ensuing statistics are not ordinary. The 
same aspect might emerge also from within the adoption 
of the (N — > oo)-perspective. We have two examples 
in mind, where the thermodynamic level, if it exists, is 
closely related to anomalous rather than ordinary sta- 
tistical mechanics. The first is the dynamical process 
recently discussed by Latora et al. ||. This is a Hamil- 
tonian system of TV rotors, each rotor being coupled to 
all the other rotors of the system with the same cou- 
pling. The authors show that in the limit TV — > oo from 
this Hamiltonian picture a process of superdiffusion un- 
der the form of a Levy walk emerges. On the other hand, 
we have in mind a second example, of Hamiltonian na- 
ture, resulting in Levy diffusion. This example is afforded 
by the class of billiards discussed in the recent book of 
Zaslavsky ||. These are systems with only one degree 
of freedom, thereby seemingly departing from the condi- 
tion TV — > oo that, according to Lebowitz jjj, is essential 
for the emergence of thermodynamics. In this paper we 
shall try to prove that the point of view advocated by 
Lebowitz Jjlf is compatible with that sustained by Za- 
slavsky |||6|,|7| in the case of fully chaotic systems, where 
the dynamical approach to statistical mechanics yields 
the same ordinary form of statistical mechanics as the 
(TV — » oo)-perspectivc, thereby making apparently use- 
less the adoption of a dynamic view. However, we shall 
argue that the role of dynamics might become important 
when the system under study is characterized by persis- 
tent memory. We make also the conjecture that a kind of 
statistical equivalence exists between extended memory 
and long-range interactions. 

For the time being, we would like to stress that in his 
very interesting book || , as well as in the original papers 
HQ, Zaslavsky focuses his attention on the anomalous 
thermodynamic nature of these generators of anomalous 
diffusion. He shows that two chaotic billiards, coupled to 
one another through a small hole in the wall separating 
one billiard from the other, result in a so strong viola- 
tion of the condition of equal distribution as to suggest 
the occurrence of a Maxwell's Demon effect. Maxwell's 
Demon was described by Maxwell in his book "Theory 



of Heat" (1871) and represented a rather simple device. 
Two chambers, A and B, are separated by a division 
with a small hole. Positioned at the hole the demon 
allows swifter molecules to pass from chamber A to B 
and slower molecules to pass from B to A. After a while 
chamber A will contain mainly swifter molecules, whereas 
chamber B will contain mainly slower molecules. This 
implies the appearance of temperature difference with- 
out expenditure of work and contradicts the Second Law 
of Thermodynamics. The version of Maxwell's Demon 
realized in practice in Refs. |^J|||f7| is equivalent to the 
case where particles with the same kinetic energy tends 
to spend more time in one of the two chambers even if 
the volume available is the same. 

According to Zaslavsky this surprising behavior is pro- 
voked by the fact that the distribution density Pii(t) de- 
viates from the thermodynamic condition of Eq. ([!]) . This 
is so because even in the phase space of seemingly chaotic 
systems there might exist stable islands. The surface of 
separation between the chaotic sea and the stable islands 
is characterized by fractal properties and these proper- 
ties make these surfaces very sticky. Consequently, in the 
long-time limit the function Pii(t) becomes proportional 
to ifj{t), denoting the distribution of waiting times at the 
border between chaotic sea and stability island. Using 
renormalization group arguments it is shown || that 



, s const 
hm = — 77T- ■ 



(5) 



where fx > 2, in accordance with the Kac theorem M. 
This means that the first moment of the waiting time 
distribution is finite. After revealing the sources of devi- 
ations from the ordinary prescriptions of thermodynam- 
ics and statistical mechanics, Zaslasky does not rule out 
the possibility that even in this case a thermodynamic 
perspective applies. He states || that all this is obliging 
us to rethink the foundation of thermodynamics. 
We see that there are several issues to discuss: 

(i) Is there a connection between the condition of 
strong chaos leading to thermodynamics in the sense of 
Eq. (|]) and the Boltzmann condition TV — > oo? Is there, 
in general, a connection between the perspective based 
on deterministic chaos of low-dimensional systems and 
the perspective based on TV — > oo? 

(ii) As far as the connection between ip(t) and the sen- 
sitivity to initial condition is concerned, we see that in 
the case where the exponential picture does not apply 
Eq.([3|) splits into the following two different possibilities: 



and 



cx £(-«) 



£(*)' 

Which is the correct property? 



(6) 



(7) 
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(iii) Is there a form of thermodynamics involved in this 
case? If it is so, which is the nature of this form of ther- 
modynamics? This would answer the fundamental ques- 
tion raised by Zaslavsky about the origin of his Maxwell's 
Demon effect §,|,§. 

(iv) Is it possible to find a general analytical expres- 
sion for ip(t), fitting the inverse power law of Eq.(||), and 
collapsing, for a critical value of a given parameter, into 
the exponential form of Eq. ([!]) with the same coefficient 

To give the reader a preliminary understanding of ques- 
tion (iv) we have to make some remarks. We shall show 
that the Manneville map Jlfj[| is a dynamic process general 
enough as to reproduce both the exponential property of 
Eq.(|l|) and the inverse power law of Eq.(||), depending 
on the value of a given control parameter z. We shall 
study this dynamic system while keeping in mind the 
non-extensive thermodynamics introduced by Tsallis in 
1988 111]] , which is now becoming more and more popular 
|^2| . From a formal point of view, within this new form 
of statistical mechanics a crucial role is played by the so 
called g-exponential. If we express the function ip(t) as a 
g-exponential we obtain 

^) = [l-(l-<7)M 1/(1 - 9) - (8) 

The parameter q is called entropic index. We see immedi- 
ately that with q > 1 the function ip(t) is an inverse power 
law. At the critical value q = 1 the function ?/>(£) would 
make an abrupt transition from the inverse power law 
Eq. (||) to the exponential form of Eq. (Q) . An important 
motivation for this paper has been to assess whether or 
not the non-extensive thermodynamics of Tsallis might 
be a satisfactory solution of question (iii) as well as of 
(iv). ^ 

This paper is organized as follows. In Section II we 
discuss the dynamical properties of the Manneville map. 
We shall discuss theoretically and numerically the wait- 
ing time distribution tp(t) in this specific case. In Sec- 
tion III we illustrate very simple arguments supporting 
the choice of Eq. (0) . Section IV is devoted to exploring 
the possibility that the Tsallis non-extensive thermody- 
namics might shed light on the thermodynamic nature of 
the Maxwell's Demon effect. Section V shows that the 
extended regime of transition from dynamics to thermo- 
dynamics corresponds to the steady process of memory 
erasure accompanying the dynamic realization of Levy 
processes. Finally Section VI is devoted to making a 
balance on the answers to questions (i)-(iv) that we are 
giving throughout this paper. 

II. THE MANNEVILLE MAP 

In this section we illustrate some key dynamic prop- 
erties of the Mannneville map pOj . The reader can find 



surprising that after mentioning the key issues of Section 
I, we decide to focus our attention on a non-Hamiltonian 
system, with only one variable. We feel therefore the 
need of justifying our choice. First of all, we want to re- 
mark that, as proved by the authors of Ref. (5), the use 
of a Hamiltonian system with infinitely many degrees of 
freedom does not rule out the possibility that the re- 
sulting dynamics, as far as a few collective variables are 
concerned, is equivalent to a random walk process that 
can be obtained dynamically also from a map ]l3| ]. For 
the emergence of low-dimensional chaos from a Hamilto- 
nian picture with a large number of degrees of freedom 
we refer to the discussion of Tennyson et al. jl4| . 

Then there is another reason why the choice of the 
Manneville map is convenient to mimic the Hamiltonian 
properties responsible for the birth of the Levy diffusion. 
According to Zaslavsky, the Levy superdiffusion in one 
of his two-dimensional systems is generated by the fact 
that the states of uniform motion last with a waiting 
time distribution that has an inverse power law nature. 
This means that any condition of uniform motion cor- 
responds to a stable island whose border with the sur- 
rounding chaotic sea is fractal. A trajectory sticks to this 
border with an inverse power law distribution of wait- 
ing times. For the whole sojourn time of the trajectory 
on this border, in the laboratory frame of reference the 
system undergoes a uniform motion in a given direction. 
From time to time the trajectory leaves the fractal region 
and after a short sojourn in the chaotic sea it sticks again 
to another fractal region. The particle sticking to a given 
fractal region moves in the laboratory frame of reference 
with constant velocity. The change from one sticking con- 
dition to another might have the effect of changing the 
motion direction, but it does not affect the velocity mod- 
ulus. Note that the velocity modulus is constant because 
the trajectory motion occurs in the constant energy sur- 
face of the phase space. A nice dynamical model sharing 
all these properties is the egg-carton model of Geisel et 
al. []l5| , where the trajectory moves of uniform motion 
along a given channel of the egg-crate landscape and at a 
given time, randomly selected from an inverse power law 
distribution, begins moving in one of the two opposite 
directions of a orthogonal channel. Also this direction of 
motion is randomly selected. Notice that randomly here 
refers to a property of deterministic chaos. 

The one-dimensional version of this generator of 
anomalous diffusion corresponds to a velocity £ with only 
two possible values, W and —W. These two states must 
have a distribution of sojourn times with an inverse power 
law. Thus, we see immediately that the one-dimensional 
version of the Hamiltonian systems studied by Zaslavsky 
becomes similar to the one-dimensional map of Ref. [[l| . 
This map has been widely used in the recent past to 
derive Levy processes |lu-0]. On the other hand, for 
the purpose of the discussion of the present paper, as 
already done in the earlier work of Ref. pfj, it is con- 
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venient to shift our attention from the bimodal maps of 
Refs. jn||l6] ^9) to the case of the Manneville map, which 
has only one laminar region. The Kolmogorov complexity 
of the Manneville map has been already analyzed years 
ago by Gaspard and Wang . The complexity of a se- 
quence equivalent to that of the bimodal maps of Refs. 
0,|6-19| was studied recently by Buiatti et at [?]. It is 



straightforward to show that the Kolmogorov complexity 
of the Manneville map is proportional to the complexity 
of the bimodal maps generating Levy processes |23] . Wc 
are interested in sporadic randomness, and for our pur- 
poses it is irrelevant whether after exiting from one lami- 
nar region and crossing the chaotic region, the trajectory 
is injected back into the original laminar region, or into 
a new laminar region equivalent to the original |23f| . If 
we make the arbitrary decision that crossing the chaotic 
region always implies a transition from one laminar re- 
gion to the other, the resulting complexity is the same 
as that of the Manneville map, with only one laminar 
region, and consequently no uncertainty associated with 
its choice. In the bimodal case, the injection back to the 
original laminar region and the jump into the other lam- 
inar region are random processes, with the same prob- 
ability. This implies a Kolmogorov complexity which is 
twice as that of the Manneville map p3| . 

In conclusion, the reason for the choice of the Man- 
neville is that this map mimics very well the condition 
of sporadic randomness that, according to Zaslavsky ||, 
is responsible for Maxwell's Demon effect. Therefore, we 
are convinced that a numerical treatments of the Hamil- 
tonian systems with sporadic randomness would lead to 
the same conclusions as those reached by us here, on the 
basis of the Manneville map. 

The Manneville map reads: 



z n (mod.l)(l < z). 



(9) 



We note that at z = 1 the Manneville map becomes 
equivalent to the Bernouilli shift map used by Zaslavsly 
to prove the fundamental result of Eq. ([[J). For z > 1 
the interval [0, 1] is divided in two regions, the laminar 
region, [0, d(z)], and the chaotic region, [<i(z),l], with 
d(z) defined by 



d(z) +d(z) z = 1 



(10) 



We review here the arguments used by Geisel and 
Thomae Jl3[ to derive an analytical expression for the 
distribution of the times of sojourn in the laminar re- 
gion. First of all we assume that the injection point Xq 
is so close to x — as to replace eq.(||) with: 

dx/dt = x z . (11) 

Thus we obtain the following time evolution: 

x(t) = [x l -* + (l-z)i\W-*l (12) 



Hence the time necessary for the trajectory to get the 
border x = d{z) is given by 



t = T(x Q ) = (- 



1 



1 



1 



z- 1 



(13) 



Note that in the special case where the initial conditon is 
so close to x — as to fulfill the condition xq << d, the 
exit time T(xq) can be satisfactorily approximated by 



T(x ) 



1 



1 



1' 



(14) 



which, as we shall see in Section III, can be used to de- 
fine the time at which we lose control of the trajectories 
departing from a region of the map very close to x = 0. 

Note that the distribution function ip(t) is related to 
the injection probability p{xq) by 



ip(t)dt = p(xo)dxo. 



(15) 



Assuming equiprobability for the injection process we 
have: 



p(x )dx 



1 dxo 



Thus we obtain: 



d{z) 1 dT 



1 dxo 



dT. 



d{zy dt 



(16) 



(17) 



By differentiating t with respect to Xq we finally arrive 
at 



^l + d*- 1 ^-!)*]-*^*- 1 ). 



(18) 



It is important to observe that the mean waiting time 
is given by 



T„ 



1 1 



(19) 



00 for 2 < z 



with T a 

We see that Eq. (|18|) implies that the region correspond- 
ing to z > 2 is characterized by a diverging first moment. 
This region is in conflict with the Kac theorem J^] and 
for this reason we do not take it into account. The region 
pertaining to the interval 1.5 < z < 2 is characterized by 
a finite second moment and a diverging second moment. 
This is the region of interest for us, since it corresponds 
to that generating Levy diffusion according to the recent 
work of Ref. |20|. The region 1 < z < 1.5 is characterized 
by a finite second moment. Of course, the ideal condi- 
tion z — 1 implies all the moments to be finite. This is 
a region which would correspond, in the perspective of 
Ref. p0[ |, to ordinary Gaussian diffusion. We note that 
in the region 1 < z < 1.5 the waiting function distri- 
bution must make a transition from the inverse power 
law behavior of Eq.dlq) to the exponential regime, where 
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the arguments of Zaslavsky leading to Eq.([y) apply. In 
fact, at z = 1, the Manneville map becomes identical to 
the Bernouilli shift map. In that case, the theoretical 
remarks of Zaslavsky yield: 



il>(t) oc exp(—tln2). 



(20) 



We are now in a position to make a preliminary ob- 
servation about question (iv) of Section I. The analytical 
expression of Eq. (|lj) can be expressed in the form of the 
g-exponential of Eq.(|h, with 



q=l + (z- l)/z. 



(21) 



However, this would not fit completely the request of 
being an exact representation of the waiting function in 
the fully chaotic case. In fact q = 1 would correctly imply 
z = 1, but this would make Eq.(|l8|) collapse into 



ip(t) oc exp(-t). 



(22) 



rather than in the exact expression of Eq. ( ]20| ) . 

This, in principle, would not rule out the possibility 
of finding a different analytical expression fulfilling the 
requirement of question (iv) . To explore more deeply is- 
sue (iv) we use a numerical treatment of the Manneville 
map of Eq. (^). For computational purposes we found 
to be more convenient to evaluate the population of the 
laminar region, M(t), rather than the waiting time dis- 
tribution ip(t). The two functions are related the one to 
the other by 



M(t) 



1 



i)(t')dt'. 



Thus, the analytical expression of Eq.(jl8|) yields 

M(t) = [i + d z -\z-i)ty 1/{z - 1) . 



(23) 



(24) 



The explicit expression of M(t) in the fully chaotic case 
(z = 1) is derived from Eq. ([!]). In fact, Zaslavsky proves 
§ that Eq.(0) becomes exact in the case of the Bernouilli 
map, with hxs — ln2. In the normalized form, Eq.(l) 
becomes: 

ip(t) = h KS exp{-h KS t). (25) 
Let us plug Eq.@ into Eq.Q. We obtain 

M(t) = exp(~tln2). (26) 



We know that Eq.(26) is exact and we expect Eq.fl24|) 
to become exact for z — » oc. In conclusion, as to the 
important function ip(t), we know its exact expression 
at z = 1 and in the asymptotic limit z — > oc. It is 
worth exploring the intermediate region with a numerical 
treatment. 

Figs. 1 and 2 illustrate the results of our numerical 
treatment. 




FIG. 1. M(t) as 3 function of time. The meaning of 

the six full lines is as follows. The lowest full line is the function 

M(t) = exp(—tln2). All the other full lines denote the long-time 

i 

inverse power law M{t) = 1/t a_1 . The dotted lines are the nu- 
merical result. All the full lines but the lowest have been shifted 
to the right to make them distinguishable from the numerical re- 
sult. The value of the parameter z, from the bottom to the top 

is: 2 = 1,1.1,1.2,1.3,1.4,1.5 




FIG. 2. M (t) as a function of time. The meaning of the four 
pair of lines is as follows. The full lines denote the function M(t) of 
Eq.(^) and the dotted lines denote the numerical results. To make 
he full lines distinguishable from the dotted lines we shifted them to 
the right by the quantity e = 0.1. In the logarithmic representation 
adopted, this is equivalent to replacing t of M(t) with texp(-e). The 
value of the parameter z from the bottom to the top changes as 
follows: 2 = 1.5,1.6.1.7,1.8,1.9 

In Fig. 1 we illustrate the result of the numerical cal- 
culation with the parameter z in the interval [1, 1.5]. We 
see that the long-time limit fits for all values of z con- 
sidered but z = 1, the theoretical prescription of the in- 
verse power law t~~^ . However, this inverse power law 
regime is reached after an extended transition regime, 
which seems to be exponential-like. The duration of this 
transition regime become more and more extended with 
z coming closer and closer to z = 1. At z = 1 this tran- 
sition regime becomes infinitely extended and coincident 
with the theoretical prediction of Eq. (|2^) . 
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In Fig. 2, devoted to studying the relaxation of the 
Manneville map with z in the interval [1.5,2], we see 
that the prediction of Eq. ( p^ ) is very accurate and tends 
to become exact with increasing the value of z. 

In conclusion, these numerical results prove that the 
q-exponential of Eq.(^) cannot satisfactorily fit the con- 
ditions set by question (iv). However, this is not yet a 
strong reason to rule out the non-extensive thermody- 
namics of Tsallis jll 12 1 . Let us make another attempt 
at establishing this non-extensive thermodynamics as the 
proper form of thermodynamics that Zaslavsky is look- 
ing for. We now set a request more flexible than the 
one earlier adopted. This is based on disregarding the 
extended regime of transition to the inverse power law. 
In other words, we now accept an analytical form of q- 
exponential which at q = 1 becomes identical to Eq.(|20|) 
while at z > 1 results in the same inverse power law as 
the Manneville map. We try with the following possibil- 
ity: 

m = [l + Q.-q)h K st] l/ll - q \ (27) 

supplemented by Eq. (pi"|) . To a first sight, this possibility 
looks very attractive. In fact, at q = 1, the prescription 
of Eq.(pl|) yields z = 1. This means the Bernouilli map, 
and, consequently B 



Kks = ln2. 



(28) 



Thus, the requirement of recovering Eq.(20) is fulfilled 



while the request of producing the correct power law at 
z > 1 is obviously satisfied due to the Eq. (pil|) . 

However, this is not physically correct. We let the 
reader know that years ago, Gaspard and Wang H] 
evaluated the Lyapunov coefficient of the Manneville 
map and found it to be given with a good accuracy by 
has = (2 — z)zln2. This means that hpcs vanishes at 
z = 2, thereby yielding a strong disagreement between 
Eq.© and Eq.(|l|). 

One might think that the departure of Eq.(^Tj) from 
Eq. (|l8|) is caused by the fact that the coefficient h q ap- 
pearing in Eq.(j^) cannot be identified with the ordinary 
Lyapunov coefficient. In fact, according to the Pesin the- 
orem [H we have 



IKS 



X 



dx p{x)ln(d<& I ' dx) 



(29) 



One might think that the adoption of Tsallis non- 
extensive thermodynamics implies that the ordinary def- 
inition of Lyapunov coefficient cannot be maintained and 
that consequently h q cannot be identified with the prop- 
erty of Eq. (p9|) , which is in fact the property discussed 
by Gaspard and Wang pl| . In Section IV B we shall 
prove that also this possibility has to be ruled out. At 
the same time we shall show that the correct answer to 
question (ii) cannot be given by Eq.(j|). The derivation 
of the correct answer to question (ii) is given in Section 
III. 



III. SENSITIVITY TO INITIAL CONDITIONS 

The sensitivity to initial conditions of the Manneville 
map has been studied in Ref. (20) . The authors of this 
paper evaluated the function £(t) defined by Eq.(||), using 
the same assumptions as those adopted in Section II to 
derive Eq. (Uq). The result is given by: 



e(t) = [i - (z - i)x*-H] 



-z/(z-l) 



(30) 



This function is, in a sense, faster than the exponential, 
since it diverges at the time Tdi V given by 



1 



{z - l)x z Q 



(31) 



while the case of exponential sensitivity would yield finite 
values of ^{t) for any finite value of time. 

It would be tempting to conclude that the Manneville 
map is characterized by a form of sensitivity to initial 
condions stronger than the usual exponential sensitivity, 
but it would be wrong. Let us explain why it is so using 
heuristic arguments. Let us divide the phase space of 
the Manneville map into two regions. The first is very 
close to x = 0, and it is characterized by the the local 
Lyapunov coefficients 



X(xo) = ln{l + zxq ) « zxq 



(32) 



The second region includes the whole chaotic region plus 
a portion of the laminar region and it is characterized by 
local Lyapunov coefficients of the order of 



A w ln(l + z). 



(33) 



Let us now define the time at which we lose control of 
a trajectory, in the case of exponential sensitivity, as a 
time proportional to 1/A, where A is the Lyapunov coef- 
ficient. We see that the time at which we lose control of 
the trajectories departing from the first region is of the 
order of the time scale Tdiv, of Eq.(|3l"|), while the time at 
which we lose control of the trajectories departing from 
the second region is of the order of unity. We note that 
for values of xq very close to x — we can fulfill the 
inequality: 



(34) 



where T av is the waiting mean time defined by Eq.(fl9|). 
Note also that for xo so small as to fit the condition of 
Eq. ( |34"| ) , the time Tdiv can be interpreted in two distinct 
ways. The former is the time at which diverges, and 
the second is the inverse of the local Lyapunov coeffi- 
cient of Eq. (|32]) . We see, in conclusion, that the com- 
plexity of the Manneville map is decided by the time 
scales l/7o<7(l + z) and T av rather than by Tdi V - This is 
the reason why the randomness detected by the method 
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of Kolmogorov complexity pi) is strongly determined by 
the influence of the chaotic region of the Manneville map. 

This is a key point worth of a comment. The function 
£(t) with the form of Eq. (p0|) cannot be used for the sta- 
tistical interpretation of the dynamic process under study 
for times larger than the exit times, and consequently 
cannot have any relevance when the time T av of Eq. ( |l9| ) 
is finite, namely for z < 2. In this case, in the time scale 
t >> T av the analytical form of the function £(t) changes 
and becomes exponential |21|. At this stage there are no 
more possible connections between tp(t) and £(t) left. In 
fact,the function ip(t) has a meaning completely different 
from that of f(t). The function ip(t) signals the proba- 
bility of observing sojourn times of any length, including 
those exceeding the time scale T av . Thus, the function 
ip(t) maintains its inverse power law nature even in the 
time scale where is exponential. 

In accordance with this remarks we find a connection 
between the function ip(t) and the function at the 
exit time T using Eq.(Q). The exit time depends on the 
initial condition of the trajectory within the laminar re- 
gion and is given by Eq. (|TJ). Let us express xq of 
Eq.(|30"|) in terms of this exit time, and let us replace it 
into Eq. ([5(j|). The resulting expression has to be inter- 
preted as the value that the function gets at the exit 
time. Its explicit expression is: 



£(T) = d^d 1 -* + {z- lJT]*/^- 1 ). 



(35) 



Comparing Eq. 
that: 



to Eq.(Uq) we reach the conclusion 



ip(T) oc 



(36) 



Thus with very simple arguments we are led to conclude 
that the generalization of Eq.(||) implies the choice of 
Eq. (^) . The meaning of this conclusion is that the value 
of ip(T) is determined by that of the function £(t) with t 
being a given exit time. This very simple argument set- 
tles question (ii) of Section 1. We have to address now 
the much more delicate question of the connection with 
the thermodynamic perspective. Before ending this sec- 
tion, we would like to notice that the condition of Eq. d36| ) 
does not apply only in the region z > 1.5 which favors 
the continuous approximation behind the analytical re- 
sults of Eqs.(|3C|) and ([l8]). It applies very well also to 
the long-time limit in the region 1 < z < 1.5, as clearly 
illustrated by Fig. 1. 



IV. THERMODYNAMIC APPROACH TO THE 
WAITING TIME DISTRIBUTION 

Here we show that the entropic arguments of the 
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would naturally yield 



conclusions of Section III. Thus the entropic arguments 
yield a wrong conclusion, which by itself does not nec- 
essarily imply that the Tsallis non-extensive thermody- 
namics is ill founded. If the Tsallis non-extensive ther- 
modynamics were well founded, on the other hand, this 
wrong conclusion would certainly imply that this form of 
thermodynamics QJl^l is not permanent. In Section V 
we shall show, on the other hand, that the Maxwell's De- 
mon effect is permanent and takes place in a later time 
scale . 



A. Entropic derivation of tp(t) 

The dynamic derivation of Levy processes from inter- 
mittent maps fl^,[l{| has established that probability for 
the particle to make a jump of length \x\ in a time t, 
II(x), is related to ip(t) by 



n(x) = ip(x/w)/w. 



(37) 



To make immediately evident to the reader the expla- 
nation of the form of Eq.(|37|) it is enough for us to let 
him/her know that the dynamic derivation of of Levy 
processes @,|9| rests on a one-dimension motion with 
a fluctuating velocity. The particle velocity fluctuates 
among the values W and — W. The particle sojourns in 
these two states with a distribution of waiting times given 
by ip(t), which has the same structure as that of Eq. (fl8|) . 
This accounts for Eq.(37). On the other hand, on the 
basis of the increasing interest for the subject of non- 
extensive entropy |12|| , and along lines similar to those 
adopted by earlier work p5|~|27t, we decide to maximize 
the Tsallis entropy to determine the form of U(x). 
Therefore we define 



s q [n(x) 



(1-9) 



(38) 



and we maximize it under the constraint of assigning a 
fixed value to the first moment of the escort distribu- 
tion with entropic index q. This means that, as done in 
Rcf . J2(| , we have to maximize the functional form F q (II) 
defined by 



F q (n) 



i-j^ujxydx 
(1-9) 



\x\n(x) q dx + a / U(x)dx. 

-oc J — oo 



As a result of the entropy maximization we get 

.4 



U(x) = 



[1 + /%-!) la^p/te- 



i) ' 



(39) 



(40) 



Eq.fl6J) rather than Eq.([1), creating a conflict with the 



It is evident that for this entropic argument to be com- 
patible with the statistical derivation of tp(t) as result- 
ing from Eq.(18|) we are forced to accept the relation of 
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Eq.©. In the region 1 < z < 1.5, as shown in Fig. 1, 
the time regime of validity of the inverse power law ful- 
filling Eq. (^l|) is confined to regions of larger and larger 
times as z — ► 1, since the size of the early time region 
with exponential-like behavior becomes more and more 
extended. 

We notice that the Tsallis information approach af- 
fords a formal justification for the birth of an inverse 
power law. However, the special value of Eq.(|2l]) is ap- 
parently established in such a way as to make the entropic 
argument compatible with the dynamic argument, with 
ad hoc arguments, rather than as a result of a theoreti- 
cal prediction. The role of the entropy of Eq. (^) would 
become much more important if it were possible to use 
it with no recourse to ad hoc arguments. In Section IV 
B we shall show that perhaps it is possible to establish a 
dynamic approach to this non-extensive thermodynam- 
ics, and consequently to a q > 1 established with no ad 
hoc assumption, with the caution, however, of interpret- 
ing this thermodynamics as temporary. We shall discuss 
also in which sense the entropic index given by Eq.(|2T 
as to be thought of as "magic" . 



B. A transient form of Kolmogorov-Sinai entropy 

The results found years ago by Gaspard and Wang 
pl[ that the ordinary Lyapunov coefficient is finite in 
the range 1 < z < 2 seem to rule out the possibility of 
applying to the present case the prediction of Eq.(|l|). In 
fact, since hxs is finite, we can use Eq. (jl|) , but this would 
lead to an exponential ip(t) in deep contrast with the fact 
that is an inverse power law instead. The breakdown of 
this important theoretical prediction is certainly due to 
the fact that the demonstration made by Zaslavsky M is 
incompatible with the existence of sporadic randomness. 

Here we plan to discuss whether this connection be- 
tween ip(t) and the KS entropy can be extended provided 
that a form of non-extensive KS entropy is used to exam- 
ine the regions of times shorter than the mean waiting 
time T av . To discuss this important possibility we must 
review the recent theoretical result of Ref . |g3 . This pa- 
per provides arguments for the non-extensive generaliza- 
tion of the Pesin theorem As shown in Ref. |22], the 
proper non-extensive generalization of the Kolmogorov- 
Sinai entropy should rest on the calculation of 



H q (N) 



i p(u ...LU N ^ 1 ) q 



(41) 



where p(u)q...u)n~i) is the probability of finding the cylin- 
der corresponding to the sequence of symbols uiq...u>n-i 
H . In the case where we set q = 1 the entropy of Eq. ( f4l| ) 
becomes the ordinary Kolmogorov-Sinai (KS) entropy. 
This entropy expression affords a rigorous way of definin- 
ing the earlier introduced concept of true entropic index 



Q. If it exists, Q is the value of the entropic index q 
making the entropy of Eq.(^l|) increase linearly in time. 

To point out the the importance of the work of Ref. 
p4j , we have to mention the difficulties with the numer- 
ical calculation of Ref. . These authors tried to eval- 
uate the dependence on N of the entropy of Eq.(^lj) to 
establish Q in the case of a text of only two symbols, with 
strong correlations. They have been forced to maintain 
their numerical study within the range of windows of size 
N = 10. The case of many more symbols would be be- 
yond the range of the current generation of computers. 
However, when the sequence of symbols is generated by 
dynamics, as in the case here under study, and the func- 
tion x) of Eq. (||) is available (for convenience, we 
keep now explicit the dependence on the initial condition 
x), it is possible to adopt the prescription of Ref. 
which writes H q (t) of Eq. Q) as 



H q (t) EE 



_ i - 51- 1 j dx P (x) q t(t, xy- q 



1 



(42) 



where the symbol t denotes time regarded as a continuous 
variable. In fact, when the condition N » 1 applies, it 
is legitimate to identify N with t. The function p(x) 
denotes the equilibrium distribution density and 5 the 
size of the partition cells: According to Ref. [^4| the phase 
space, a one-dimensional interval, has been divided into 
M = 1/(5 cells of equal size, with M >> 1. 

If the property of Eq. ( |30| ) held true with no dependence 
on the initial condition, namely in the form 



at) = [i-(z-i)\ z t]-*^- i \ 



(43) 



one would reach the conclusion that the entropy of 
Eq. ( [l2"| ) increase linearly in time if we adopt the entropic 
index of Eq. (pij ) . At the same time the theoretical pro- 
posal of Eq. (||) would turn out to be correct. This is the 
reason why we call this special value of q 'magic" and we 
denote it with the capital letter Q thus writing 



Q = l + (z-l)/z. 



(44) 



What about the case under study, where the local Lya- 
punov coefficient depends on xl A possible way out 
would be given by the following procedure. We assign to 
the entropic index q of Eq.jJL?]) the value prescribed by 
Eq. . This has the nice effect of making the entropy 
H q {t) linear dependent on time. Unfortunately the distri- 
bution function p{x) raised to the power of 1 + (z — l)/z, 
results in a divergent contribution as it can be easily 
checked noticing [^l| that p(x) oc l/x z ^ 1 . We notice 
that it would have been very appealing to establish the 
"thermodynamic" foundation of ip(t) on the relation 



m =< (a-t)) 1 



> 



i/(i-Q) 



(45) 



with the mean value < A(x) >q defined by 
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< A(x) > Q = / dxp(x) Q A(x) 



(46) 



This possibility is ruled out by the divergency stemming 
from p(x)Q. 

We have to point out another important property sug- 
gested by Eq.(H|). We see that this formula shows that 
the only possible way of making H q {t) independent of 
the cell size 8 is to set q = 1. In Ref. |28| it has been 
pointed out that in the multifractal case it is convenient 
to carry out the average over the distribution of power 
indices rather than on p(x). This has the nice effect of 
making H q (t) independent of the cell size. In the spe- 
cific case where the invariant distribution exists and is 
not multifractal, as in the case here under study, we are 
compelled to use q = 1. This is line with the fact that 
when the invariant distribution of the Manneville map 
exists, namely if z < 2, the average Lyapunov coefficient 
is finite, the Pesin theorem work and the ordinary KS 
entropy can be defined . 

Therefore we are left with the possibility of applying 
the prescription of Eq.([42|), with the associated average 
of Eq.(f46|), to the case where p{x) is not the invariant 
measure, but it is a given out of equilibrium distribution. 
This would correspond to giving up the assumption that 
the Tsallis non-extensive thermodynamics has an equi- 
librium significance, in line with the numerical results of 
Rcf. [§9). 

On the other hand, it is evident, on the basis of the 
conclusions of section III, that the function £,(t,x) does 
not keep its power-law nature forever, and that in the 
time scale t >> T av , with T av given by Eq.(|l9|), it is ex- 
pected to recover an exponential form. Thus, are back to 
the point that this non-extensive form of thermodynam- 
ics must be temporary. 

We would be tempted to conclude this section with 
the following answer to question (iii) of Section I. The 
Maxwell's Demon thermodynamics is the non-extensive 
Tsallis thermodynamics, and this is a form of out of 
equilibrium thermodynamics, implying that also the 
Maxwell's Demon effect is temporary. 

This conclusion, to a first sight, seems to be reasonable. 
However, it yields a conflict with the result of the earlier 
work of Ref. E^] . The authors of Ref. [^2| studied the en- 
tropy of Eq.(|4l|) with a symbolic method, and, as earlier 
said, they have been forced to use short-time windows. 
The adoption of short-time windows is imposed by the 
fact they did not make any direct recourse to dynamics as 
in the calculation yielding Eq. ([f3|) . The case studied in 
Ref. [^2| refers to a dynamic system with the same com- 
plexity as the Manneville map here under study. Proba- 
bly the main difference with the analysis here illustrated 
is that the method adopted by the authors of Ref. p^ ] 
is equivalent to studying the Geisel map rather than the 
Manneville map, with a repartition of the space into only 
two cells. The conflict between the analysis here made 



and the result of the work of Ref. |2^] is probably due to 
the fact that here we adopt a much finer repartition of 
the phase space. 

This kind of disagreement between two different theo- 
retical predictions about the entropic index of the same 
dynamic process is not yet a reason for the rebuttal of the 
non-extensive thermodynamics as the explanation of the 
Maxwell's Demon effect. In Section V we shall see that 
the Maxwell's Demon effect seems to be caused by the 
emergence in chamber A 7 made more densely populated 
by the Demon, of a diffusion equation with a "diffusion 
coefficient" weaker than that of the diffusion equation 
emerging in chamber B. This dynamic effect is explained 
using the dynamic derivation of Levy statistics, which im- 
plies memory erasure and consequently emerge in a time 
scale where the temporary non-extensive thermodynamic 
regime, if it ever exists, is over. 



V. LEVY PROCESSES AND MEMORY ERASURE 

Lebowitz |ij states that the Universe time evolution is 
characterized by the transition from a condition of lower 
to a condition of larger entropy without implying a de- 
parture from microscopic reversibility. Let us focus our 
attention on the diffusion process: 



dx , . 



(47) 



Let us assume that the stochastic velocity v(t) fluctuates 
among the two values W and —W, and that the distri- 
bution of times of sojourn in one of these two velocity 
states, the same for the two states, is given by Eq.(Q. 
It is shown |30| that the distribution density p(x, t) obeys 
the following equation of motion 



g- t p(x,t) =<v 2 > J <£ v (t')—p(x,t-t')dt', 



(48) 



where < v 2 > is the equilibrium mean quadratic value of 
v, and <fr v (t) is the normalized correlation function of the 
variable v. For details, and especially to understand why 
this equation in the long-time limit yields the process of 
Levy diffusion, the interested reader should consult Ref. 
fj. 

Here we limit ourselves to noticing a few aspects. First 
of all, the distribution density p(x, t) can be interpreted 
as a probability concerning a very large number of par- 
ticles imbedded in the same space phase, for instance, 
the phase space of the billiards of Zaslavsky Thus, 
in a sense, the condition N — » oo is fulfilled. This has 
to do with question (i) of Section I. If condition N — > oo 
is essential for the birth of statistical mechanics, what is 
the role of mixing, and that of ensuing ergodicity? 

We try to answer this question by making a second, 
very important, remark. This is that Eq.(fl8|) is exact 
under the following conditions: 
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(a) The variable v(t) is dichotomus, with only two pos- 
sible values, either W or —W. 

(b) All the particles are initially located in the same 
position, x = 0. 

We note that condition (a) is that mentioned in Sec- 
tion IV A to make easier for the reader to understand 
the meaning of the key relation of Eq.(j37|), leading to 
the entropic derivation of the waiting time distribution 
ij)(t). Both conditions are relevant for our main purpose 
of comparing the point of view of Zaslavsly fl to that 
of Lebowitz [Q. The latter property certainly realizes 
an initial condition whose probability is extremely low, 
thereby making impossible the regression to the initial 
condition. In the ideal case of a time inversion, this would 
be possible, in the irrealistic case of no perturbation of 
any kind, ranging from the round-off errors to the en- 
vironmental fluctuations. This is again in line with the 
observations of Lebowitz on the birth of statistical me- 
chanics. 

The former property, however, serves in our opin- 
ion a purpose which apparently conflicts with Lebowitz, 
namely, that of stressing the importance of mixing, as 
a key property leading to an invariant distribution. For 
the reader to appreciate this aspect it is necessary that 
he/she goes through the derivation of Eq.(|48|) JT^j. The 
projection method adopted in Ref. jl9| to derive Eq.ji^) 
rests on a projection operator, which, in turns, implies 
the existence of an invariant distribution for the variable 
v. In the specific case of deterministic chaos with mixing, 
the prototype of which is given in fact by the Bernouilli 
shift map (as we have shown, the Manneville map with 
z = 1), an invariant distribution exists and is quickly 
reached from any off-equilibrium condition. 

The waiting function ip(t) is proportional to the sec- 
ond derivative of the correlation function & v (t). This 
means that the departure from the exponential condition 
of Eq. (Q) has dramatic effects on the process of transition 
to statistical mechanics, as birth of a diffusion process. In 
fact, in the region 1.5 < z < 2 it is still possible to fulfill 
the condition necessary for the existence of the invariant 
distribution according to the Kac theorem [||. However, 
in this specific case the correlation function satisfies 
the asymptotic property 

§ v (t) = const /t p , (49) 

with 

= _L- - 2. (50) 
z — 1 

This means that in the interval 1.5 < z < 2 the index 
(3 runs between (3 = 1, at the border with the basin of 
Gaussian attraction, and (3 — 0, corresponding to the 
critical condition of no relaxation of the velocity correla- 
tion function. 

In conclusion the condition of strong chaos, with 
strongly mixing properties, has the effect of producing 



ordinary Brownian diffusion. It is not clear to us what 
would be the consequence of a dynamic condition with no 
form of mixing. We think that in principle it is possible 
to produce a fast decay of the correlation function (i) 
as a result of a mere superposition of infinitely many 
normal modes (again in action the condition N — » oo). 
However, we are not aware of any treatment proving that 
a natural equilibrium distribution is reached also in that 
case, if no arbitrary statistical assumptions are made. 

We have the impression that in the case of ordinary 
statistics it might be difficult to support the dynamic 
view against the (N — > oo)-perspective, or the latter 
against the former, without using very subtle arguments 
such as the request of neither explicit nor implicit statis- 
tical assumptions, which, in turn, quite probably would 
be questioned. For all practical purposes, the two views 
are expected to provide equivalent results. For instance, 
in the dynamic perspective adopted by Zaslavsly the 
Poincare recurrences are frequent and those with very 
short time duration are more probable than those of very 
large time duration (see Eq.([j])). This seems to conflict 
with Boltzmann's idea that the Poincare recurrences of 
systems with very large number of freedoms become ex- 
ceedingly large, so as to exceed the possibility of any 
direct observation. However, if we make the assump- 
tion that the experimental observation is made on a large 
number of non-interacting particles, all of them moving in 
one of the billiards studied by Zaslavsky PUffl 1 we reach 
a different conclusion. The phase space under study is 
not more that of the two-dimensional billiard: It becomes 
a 2N-dimension phase space. The distribution density 
driven by Eq.([48|) becomes a 2iV-dimensional trajectory, 
whose return to the initial condition takes place, in ac- 
cordance with Boltzmann's view, virtually after an in- 
finitely long time. In Zaslavsky's picture the process cor- 
responding to the 2iV-dimension trajectory leaving the 
initial condition and returning to it after an extremely 
long time would correspond to the following scenario. 
Let us imagine, for instance, that the initial distribution 
density is different from zero only in a small fraction of 
the whole 2-dimension phase space, where it is assumed 
to be constant. As a result of mixing, this initial dis- 
tribution, throughout its time evolution, would undergo 
a fragmentation process that, within a kind of coarse- 
grained perspective, would correspond to the volume of 
this initial condition growing till to becoming identical to 
the volume of the whole phase space. The Poincare re- 
currences of the Boltzmann-Lebowitz perspective would 
correspond to this dilated distribution shrinking to the 
initial condition, a process that is beyond the range of the 
experimental observation as the return of the Boltzmann- 
Lebowitz trajectory. 

We think that the relevance of dynamics for statisti- 
cal mechanics cannot be easily ruled out in the case of 
processes with long-time memory as those here under dis- 
cussion. Quite on the contrary, here we show that the dy- 
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namic approach to Levy diffusion in the form discussed in 
Ref. (3^] can account for the Maxwell's Demon effect. In 
this moment, it is not clear to us how this interesting ef- 
fect might be explained using the (N — > oo)-perspective. 
The authors of Ref. Jpfl] prove that in a time scale much 
larger than T av Eq. (|48|) as an effect of the memory erasure 
provoked by sporadic randomness becomes equivalent to 
an equation of motion that in the Fourier representation 
reads 



^p(k,t) = -b\kf +1 p(k,t), 



with 



b = T(-(3) wP +1 rip 







(0 s 

cos(--k). 



(51) 



(52) 



Note that p(k,t) is the Fourier transform of p(x,t) and 
that this is the Fourier representation of the well known 
process of Levy diffusion. The element of interest of the 
result of the dynamic approach of Ref. |30| is that the 
coefficient b of Eq. ( |52] ) plays the role of a diffusion coeffi- 
cient. The reader can easily realize that the intensity of 
this coefficient can be different in the two chambers of the 
experiment on the Maxwell's Demon effect [|||J^] even if 
the velocity W is the same. The possibility of realizing 
two different values for b using the same kinetic energy 
becomes still wider if we refer to a generator of fluctu- 
ations with an inverse power law distribution of waiting 
times more general than the Manneville map |p3"lj . It is 
evident that two different diffusion equations in the two 
chambers can produce a breakdown of the condition of 
equal population. The trajectory, in fact, will tend to 
spend more time in the regions with a smaller "diffusion 
coefficient" . 

In conclusion, according to this heuristic interpreta- 
tion, the Maxwell's Demon Effect is permanent and it 
can show up after the long time process of transition from 
dynamics to diffusion. This rules out the non-extensive 
thermodynamics of Tsallis as a possible theoretical inter- 
pretation of this effect. In fact, in Section IV B we have 
seen that if the non-extensive approach applies, it does 
as a form of temporary thermodynamics. 

We can also rule out the non-extensive nature of en- 
tropy as the basic property behind the Maxwell's Demon 
Effect. Let us see why. We have two chambers, A and B, 
separated by a wall, with a little hole. A can be a Sinai 
billiard and B a Cassini billiard |7j . The physical entropy 
of this regime, S, expressed in terms of the conventional 
Gibbs prescription, might not be extensive. The entropic 
indicator advocated by Tsallis |lj] is not additive, and it 
violates the property 



S(A + B) = S(A) + S(B) 



(53) 



even when the two subsystems A and B are totally un- 
corrected the one from the other. In the case of the 



dynamic approach to Levy diffusion the breakdown of 
the additivity condition might be dictated by memory 
effects. To explain why it is so with intuitive arguments, 
we can refer ourselves to the dynamical experiment itself 
used by Zaslavsly and Edelman [Q to reveal the existence 
of Maxwell's Demon effect. A particle that moves from A 
into B might still have memory of its initial condition in 
A, if sporadic randomness did not act for a sufficiently ex- 
tended time. This means that the transition process that 
in the long-time scale will produce Levy statistics is not 
extensive in nature. However, we expect that when Levy 
statistics are finally established, the usual additive condi- 
tion is restored, as a consequence of the memory erasure 
process. Using the earlier arguments we can imagine the 
possibility that the two chambers might correspond to 
two different difusion equations, even if the velocity in- 
tensity of the particle does not change with moving from 
the one to the other chamber, thereby resulting in the 
Maxwell's Demon effect after the process of memory era- 
sure. 



VI. CONCLUDING REMARKS 

Let us summarize our conclusion about the questions 
raised in Section I. 

Question (i). Our arguments concerning question (i) 
are not compelling and rigorous. However, we hope that 
they might make plausible the following statement: Both 
deterministic chaos, whose role is well understood only 
in the case of low-dimension systems, and the condition 
N — > 00 give rise to the birth of statistical mechanics. 
We are convinced that the joint use of these two per- 
spectives will contribute a deeper understanding of the 
dynamic origin of statistical mechanics. It is interesting 
to notice that the exponential of Eq. (|l|) , implying a fre- 
quent return to the initial condition, does not conflict 
with the irreversible character of thermodynamics and 
statistical mechanics if we perceive the Gibbs ensemble 
method as being equivalent to the (N — > oo)-perspective. 

Question(ii). We found the proposal of Eq.(H), valid 
at the exit time of a trajectory from the laminar region, 
to be the correct connection between the sensitivity to 
initial conditions and the distribution of waiting times 
in the laminr region. This means that in the case of the 
dynamical systems studied by Zaslavsky, the distribution 
of the Poincare recurrence times can be related to the KS 
entropy only in the absence of intermittency. When an 
intermittent process is present, the long-time form of the 
distribution Pn(t) changes from the exponential form of 
Eq.(^) to an inverse power law form. In the region 1.5 < 
z < 2, corresponding to the emergence of Levy processes, 
the KS entropy is still finite and the adoption of Eq.(Q) 
would conflict with the observed inverse power law nature 
of this distribution. The non-extensive thermodynamics 
of Tsallis O] might help, provided that it is viewed as 
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temporary and provided that the conflict between the 
prediction about whether Q > 1 or Q < 1 is settled. 

Question (iii). The dynamic approach to Levy diffu- 
sion of Ref. |5(| proves that the Maxwell's Demon effect is 
permanent and it takes place in a long-time scale. This 
means that it is not necessary to use memory proper- 
ties to explain it, and consequently, not even the non- 
extensive nature of entropy. 

Question (iv). The condition of ordinary statistical 
mechanics does not imply an abrupt transition when we 
set z = 1. However, for even infinitesimally small de- 
viations from z = 1 in the long time limit the waiting 
time distribution ip(t) appears to be characterized by an 
inverse power law form. This does not have strong con- 
sequences, though, because in the region z < 1.5 the 
second moment of the waiting time distribution is finite, 
and consequently the system is attracted by the Gauss 
basin of attraction. Nevertheless, the waiting time dis- 
tribution in the region where z is very close to 1, is an 
exponential at short time and an inverse power law at 
long time, a property which makes it incompatible with 
the structure generated by the method of entropy maxi- 
mization resting on the non-extensive indicator of Tsallis 

©■ 

It seems to us that the long-standing problem of a 
thermodynamic approach to Levy processes |3l| is not 
yet solved, and that the Maxwell's Demon effect of Rcfs. 
HHJl]] can be regarded as a paradigmatic case challeng- 
ing the Boltzmann perspective advocated by Lebowitz 
jl| . We plan to do further research work along these di- 
rections. 
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